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Abstract 

A solution for the neutrino mass and mixing pattern is proposed which is compatible 
with all available experimental data on neutrino oscillations. This solution involves 
Majorana neutrinos of the pseudo-Dirac type, i.e. mMajorana *C wiDirac- The s °l ar 
and atmospheric neutrino observations are mainly explained as v e — and — 
oscillations, where S indicates the sterile ( "righthanded" ) partner of each neutrino 
generation, while the LSND result is interpreted in terms of standard v^—v e oscillations. 
The resulting constraints on — v T and v T — oscillations are also discussed. This 
solution leaves room for a hierarchical mass and mixing scheme with a u T mass in 
the few eV range, as favoured by some dark matter scenarios. The apparent conflict 
with standard Big Bang nucleosynthesis is addressed and the implications for current 
and future experiments are discussed. It is argued that both short and long baseline 
accelerator neutrino experiments are needed in order to decide between this solution 
and other oscillation scenarios. 
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1 Introduction 



Currently there is a wealth of unexplained phenomena in neutrino physics which can be in- 
terpreted as indications for the existence of neutrino oscillations. The solar neutrino problem 
U, lacking a satisfactory astrophysical solution, is commonly attributed to the disappear- 
ance of electron (anti)neutrinos into some other neutrino type. Similarly, the atmospheric 
neutrino anomaly is often interpreted as an indication of neutrino oscillations involving 
Up disappearance, with or without partial reappearance as v e . The LSND experiment [[J 
claims direct evidence for — v e oscillations in a region which is partially unconstrained 
by other experiments. Furthermore, neutrinos are a prime candidate for a partial solution 
to the missing dark matter problem || if at least one mass eigenstate lies in the eV range. 
These current indications are summarized in Fig. [[]. 

Within the framework of 3 generation neutrino mixing, all recent solutions ||||@ to 
this neutrino puzzle either discard some of the experimental evidence or have to accept a 
bad fit to part of the data. Other solutions || involve the existence of one additional light 
sterile neutrino, yielding at least two extra degrees of freedom (mass and mixing angle (s)). 
Solutions involving more than one such neutrino P| ||10|| have often been discarded due to 
bounds on the number of neutrino types suggested by Big Bang nucleosynthesis [0 Ell- 
in this paper a solution is presented based on the assumption that neutrinos are of the 



pseudo-Dirac type, i.e. mMajorana -C mDirac [ 13 1 . This involves the splitting of each neutrino 



generation into two almost degenerate mass eigenstates, leading to neutrino-antineutrino 
oscillations with maximal mixing, similar to oscillations in the K° — K° system. Since oscil- 
lations do not change the handedness, left-handed neutrinos (right handed antineutrinos) will 
transform into left-handed antineutrinos (right handed neutrinos) which do not participate 
in weak interactions, and therefore appear sterile. 



In contrast to earlier schemes of this kind pOfl , the neutrino magnetic moments are 
assumed to be 0, therefore avoiding v L <-> u R transitions. Special emphasis is placed on the 
case of exactly maximal mixing, which minimizes matter effects. Furthermore, the scheme is 
extended to include 3 generation flavour mixing, and applied to the most recent experimental 
data. 

With only one parameter (a mass splitting parameter chosen to be universal for sim- 
plicity) in addition to the standard 3 generation mixing, all experimental indications for 
neutrino oscillations can be fully explained in the context of a hierarchical mass and mixing 
pattern, analogous to the quark sector. In addition, this solution suggests a substantial neu- 
trino contribution to the missing dark matter. A mechanism through which the Big Bang 
nucleosynthesis bounds mentioned earlier can be evaded is discussed in the appendix. 



2 The pseudo-Dirac neutrino formalism 



The most general neutrino mass term (Dirac-Majorana mass term) that can be added to 
the standard model Lagrangian leads, for one neutrino generation, to the (Majorana) mass 
eigenstates JT3]] JTl] 



mi ii 



(m L + m R ) ± J (m L - m R ) 2 + Am 



D 



1 



where mo is the Dirac mass, and m^,m^ are the left and right handed Majorana masses, 
and to a mixing angle 

tan 29 = 2m£>/ {tur — ttil) 

between these two states. 

The standard Dirac case is obtained by setting mi = m^ = 0, while rri£> = leads to 
pure Majorana neutrinos. The limit tjir ^> mp, m^ — (9 — m^jmn <C 1) yields the 
mass eigenstates mi = m 2 D /mn, mn ~ mn which corresponds to the well-known see-saw 



model where is assumed to be of the order of the quark and charged lepton masses, 
and m,R is a large mass of the order of some unification scale. This forces the mass of the 
"left-handed" neutrino (mi) to be small, and makes the "right-handed" neutrino (mn) very 
heavy, therefore decoupling it from interactions in the accessible energy range. 

Finally, the limit ^> m^^mn (9 ~ 45°) yields almost degenerate mass eigenstates 
mi n — mo, with close to maximal mixing arising in a "natural" way. This case, which is 
known as the pseudo-Dirac case, will be pursued further here. Taking m^ = mn = 5m for 
simplicity one obtains 

mi, ii — m D i Sm, 9 = 45°. 

In this case, the mass eigenstates can be written in terms of the neutrino/antineutrino 
eigenstates of the weak interactions as 

Wi >= -j=(W > +\z >), Wn >= -j=(W > -\* >), 

a pattern which is similar to the one in the K° — K° meson system. As in the K° case, 
this yields maximal oscillations between particle and antiparticle states if pure neutrinos or 
antineutrinos are produced. Since, up to a small mass effect, neutrinos are emitted left- 
handed and oscillations cannot change the handedness (angular momentum conservation) 
they will oscillate into left-handed antineutrinos which appear sterile. The same is true for 
right-handed antineutrinos oscillating into right-handed neutrinos. Thus phenomenologically 
one obtains active-sterile neutrino oscillations with 

Am 2 = (mD + 5m) 2 — (mn — 5m) 2 = Am^Sm 

and maximal mixing (sm 2 29 = 1). 

This can easily be generalized to 3 generations. Neglecting CP violation, the resulting 
mixing matrix is a real 6x6 matrix in principle. Here, we simplify it with the hypothesis that 
the oscillations between different generations and the active-sterile oscillations decouple, thus 
yielding the standard 3 generation mixing matrix U a i (a = e, /i, r, i = 1,2,3), equivalent to 
the CKM matrix [0, plus maximal active-sterile mixing for each generation 



\v a >= U ai -^=(\vij > +\u itII >), \v a >= U ai -^=(\vij > -\vi,n >)■ 

The above discussion remains valid to good approximation for m^ ^ mn as long as both 
terms are very small, and in particular for = 0. Since mn mass terms can be introduced 
ad hoc as bare mass terms without destroying gauge invariance, no attempt is made at 
this stage to clarify the origin of the Majorana masses. Phenomenologically, two cases are 
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considered: (1) same Srrii = e x M for each generation, where M = 1 eV is an arbitrary 
fixed reference mass, and (2) mass dependent splitting Srrii = e x mm. The mass pattern 
generated by these scenarios is therefore given by rani — mx,m2,ms for the 3 generations, 
and by the single additional parameter £ < 1 fixing the mass fine splittings. 



3 The neutrino oscillation solution 

To find a simple general solution to the neutrino oscillation puzzle described in the intro- 
duction the following further assumptions are made, in analogy to the quark sector: 

• The mass structure is hierarchical, i.e. ni\ <C 777.2 *C rriz- 

• Mixing between adjacent generations is small (#12 and #23) and mixing between far 
generations is negligible (#13 ~ 0). Flavour eigenstates can therefore be approximately 
identified with mass eigenstates (1 ~ e, 2 ~ 3 ~ r). 

Including e this yields a total of six free parameters (toj(z = 1, 2, 3), e, 9i 2 , #23)- 
The Am 2 values relevant for oscillations are Am 2 2 , Am 2 ^ for conventional mixing, and 
Am|/ JZ j = irriiMe for case (1) (Am 2 -/^, = Amfe for case (2)) for active-sterile oscillations. 
In the simplified case of only two contributing flavour and mass eigenstates, the probability 
of a neutrino of type a to oscillate into type f3 is given by the well known formula 

P a/3 = sin 2 26 aP sin 2 (1.27 Am 2 a/3 L/E) 

where L is the distance travelled in km, E is the neutrino energy in GeV, and Am^ is 
expressed in eV 2 . 

In the following, a complete solution for the six free parameters is derived based on all 
known experimental evidence (table [I]). A graphical representation of this solution is shown 
in Fig. |T[ 



3.1 LSND 



The LSND result || is interpreted as direct — v e oscillations. The large Am 2 region 
(Am 2 > 10 eV 2 ) is excluded by NOMAD |7|, CCFR @, BNL 734 @, and BNL E766 



0[ . The latter constrains also the island around Am 2 ~ 6 eV 2 advocated by Caldwell |2T| . 
The remaining region (Am 2 < 2 eV 2 ) is fully consistent with, and even slightly favoured by, 
the recent KARMEN1 limits [fj^] partially based on a small (~ la) positive result. At large 
sin 2 26 l the result is constrained by reactor experiments, including Gosgen [23J, Krasnoyarsk 
2l, Chooz 



and Bugey fl26| . Note that the Bugey result, which yields the most stringent 
limit, is also based on a small (up to 1.5 a) positive effect | 26| , fully consistent with the 
limits from the other reactor experiments. Inspired by this effect we choose a reference 
solution Am 2 2 ~ 0.5 eV 2 , sin 2 2# 12 ~ 2 x 10 -2 . However, no significant lower constraint on 
sin 2 2# can be derived from the Bugey result, allowing sin 2 2# 12 to be as small as 2 x 10~ 3 for 
Am 2 2 ~ 2 eV 2 . 
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Figure 1: Current indications for neutrino oscillations in the sin 2 26 l vs. Am 2 plane, as- 
suming dominance of two flavour oscillations in each case. The figure shows the 90% c.l. 

— > v e region favoured by LSND (double hatched), including constraints from BNL 
776 and Bugey, the Superkamiokande disappearance solution to the atmospheric neutrino 
problem (hatched), the 3 mutually exclusive 95% c.l. solutions to the solar neutrino problem 
(light shaded), and the 90% c.l. limit on — > v T (black line). Also shown is the parameter 



space of the proposed solutions for vx — > Vx (bars at sin 29 = 1) 



v.. 



v T (dark shaded) 



and — * v e (double hatched, identical to LSND region). The reference solutions described 
in the text are also indicated (black circles). 
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3.2 Solar neutrinos 



The solar neutrino problem is interpreted as v e — z/f oscillations with sin 2 2# e s = 1. This 
yields an average expected suppression factor of 0.5, slightly reduced further by a 0(1%) 
correction from v e — oscillations. Energy dependent oscillations have to be introduced 
in order to describe the observed suppression factors |27j with respect to the standard solar 
model p8| of 0.51 ± 0.05 for the Gallium experiments |2 
experiment pi, and 0.39 ± 0.03 for the water Cerenkov detectors fl3Tf |32|| 



0.27 ± 0.02 for the Homestake 
Since the MSW 



solution [|27|]f33|] is irrelevant for maximal mixing, the vacuum solution [ 27||[p4 is the only 



possible one. This yields Amj /7/1 ~ 5 x 10 - 10 u eV 



3.3 Atmospheric neutrinos 



The atmospheric neutrino anomaly is interpreted as — oscillations with sin 2 26 IJi s — 1- 
A small effect from — v e and — v T oscillations should be added. The Am 2 for v e disap- 



pearance (Amj/jf x ) is too small to have any influence, so electron-like events are expected 
to be slightly enhanced by — v e oscillations without L/E dependence. On the other hand, 
muon-like events are depleted mainly by z/ M — oscillations, with a small L/E independent 
suppression from — u e , and another small but potentially complicated contribution from 
z/ M — v T followed by v T — z/f . Using the LSND reference solution quoted above, an L/E 
independent suppression would imply a value of R = fl ^ measured Q f 0.49 or lower (if z/„ — v T is 

AV ^expected 

significant). The observed average suppression of only about 0.63 is therefore an indica- 
tion for non- maximal suppression due to dependence on L/E. Indeed such a dependence 
is observed indirectly through a dependence on the zenith angle by Kamiokande |35| and 
Super kamiokande f36|. Reinterpretation of the Superkamiokande result in terms of maximal 
— oscillations suggests a value of Am 2 ^ //2 in the range 3 x 10 -4 — 10~ 2 eV 2 |J7J. The 



detailed consideration of matter effects can potentially induce a slight upward shift of this 
range [38]. Such low Am 2 values can also explain why the suppression is stronger for the 
water Cerenkov detectors |3S| [^| than for some of the iron detectors ]39| which have a higher 



energy threshold. Experiments looking for upward going muons yield inconclusive results so 



far [40] [41], but are not inconsistent with the z/ M — interpretation. 



3.4 Dark matter 

The potential contribution of neutrinos to hot dark matter is given by Vt v = Yjm v j (92h 2 eV) 
f42|1 , where f2„ is the neutrino energy density parameter expressed as the fraction of the crit- 
ical density Q = 1 needed to close the universe, h is the dimensionless renormalized Hubble 



constant (0.5 < h < 0.85 [f4~3]) and Sm y is the sum over the masses of all neutrino types. For 
pseudo-Dirac neutrinos each neutrino flavour enters this sum twice if its sterile counterpart 
reaches thermal equilibrium in the early universe, and at least once otherwise. In the hierar- 
chical mass pattern advocated here this contribution arises mainly from r neutrinos ] 44| . No 



direct evidence for — v T oscillations exists so far in the context of the pseudo-Dirac solu- 
tion. Constraints are imposed by the exclusion regions obtained by accelerator experiments 
r7|[4li]pH][i7|. Two typical solutions close to the range accessible by current experiments 
and relevant to the dark matter problem could be (a) Am^ = 100 eV 2 (777,3 — 10 eV) 
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with sin 2 2# 23 = 10~ 3 , and (b) Am 2 3 = 3 eV 2 (m 3 ~ 2 eV) with sin 2 2# 23 = 5 x 10" 2 . 
Solution (a), with its very small mixing, would have virtually no effect on any of the current 
measurements with positive indications for oscillations, but yields a large hot dark matter 
component (Q u ~ 0.4 for h ~ 0.7) which is disfavoured by most current dark matter scenar- 
ios 0] . The associated Am 2 3 could however be lowered to any value larger than 1 eV 2 (mass 
hierarchy limit), corresponding to Q v > 0.03, without significantly changing the conclusions. 
Alternatively, it could be raised up to the extreme limit from pure hot dark matter mod- 
els of about (30 eV) 2 | j44}| . Solution (b) would yield a small but significant contribution of 
about 5% to the atmospheric neutrino suppression, with only a moderate hot dark matter 
contribution {VL U ~ 0.1 for h ~ 0.7). 

3.5 The mass fine splitting 

The Am 2 for v T —v^. oscillations is essentially unconstrained by current data. Since ni2 ~ 1 eV 
(LSND), the atmospheric neutrino solution yields e ~ 10~ 4 — 2 x 10~ 3 for both case (1) and 
case (2). Folding this with the Am 2 3 parameter space spanned by — v T solutions (a) 
and (b), Am^ JJ3 could be anywhere in the range ~ 10 -3 — 1 eV 2 . In combination with 
v n — v r oscillations this could yield a small contribution to the angular effect observed in the 
atmospheric neutrino anomaly. 

Finally, the above e value applied to the solar neutrino result predicts mi in the range 
10~ 9 — 10~ 6 eV for case (1), and 10~ 5 — 10 -3 eV for case (2). This mass and mixing pattern 
is in agreement with the nonobservation of neutrinoless double (3 decay p8| . 

Table 1: Parameters of the pseudo-Dirac neutrino solution for e = 10 -4 — 2 x 10~ 3 (cases (1) 
and (2) of section [3~5| combined), m refers to the Dirac mass of each neutrino generation and 
Am 2 to either the mass fine splitting inside each generation, or the mass difference squared 
between generations where appropriate. sin 2 2# refers to the associated mixing angle. Only 
the entries in the m and sin 2 2# columns are independent parameters, Am 2 being fixed by m 
and e. 



generation 


m (eV) 


Am 2 (eV 2 ) 


sin 2 29 


1 ("e") 


10~ 9 - 10~ 3 


5 x 10~ 12 - 10" 10 


1 


2 (V) 


0.5-1.7 


3 x 10" 4 - 10~ 2 


1 


3 (V) 


1-30 


4 x 10~ 4 - 2 


1 


1^2 




0.3-3 


0.002 -0.035 


2^3 




1-30 
or 30 - 900 


< 0.07 
< 0.0035 


1^3 




1-900 


negligible 



4 Alternative solutions 

Only two of the critical experimental inputs are so far not cross-checked independently 
by at least one other experiment: The LSND result [Q, and the Homestake result |5[J 
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for solar neutrinos. The former fixes the — v e oscillation parameters and, if disproven, 
would essentially remove any lower constraint on Am 2 2 . Comparison of the latter to the 
other solar neutrino experiments yields the main evidence for the energy- dependence of the 
solar neutrino deficit. A 30% increase of the Homestake neutrino flux measurement (not 
compatible with the quoted error), and a reduction of the predicted B s solar neutrino flux 
by about 20% (which is inside the currently discussed error margin ||49|| ) would make all solar 
neutrino measurements compatible with a flat suppression factor of about 0.5 (see also 00). 
This would remove the constraint on Am^ n 1 which, ignoring Big Bang nucleosynthesis 
implications, could then be as large as 10~ 4 eV 2 (bounded by its potential effect on the 
atmospheric neutrino anomaly). In order to fix the parameters of the pseudo-Dirac model 
the verification of LSND by KARMEN2 |22| and later by other experiments |K| is therefore 
crucial, and a verification of the energy dependence of the solar neutrino suppression, e.g. by 



Borexino ||5l| , is highly desirable. Note that the azimuthal dependence of the atmospheric 
neutrino suppression, which is observed by only one type of experiment, is not an a priori 
independent experimental input in this context, but rather a prediction of the model based 
on the observed non-maximal suppression. 



5 Predictions for current and future experiments 

The solution of the neutrino oscillation puzzle proposed in this paper can be checked by 
several current and future experiments. 

Since the solar neutrino problem is solved via v e disappearance into sterile neutrinos, 



SNO [52] should not observe an enhancement of the NC/CC ratio as it would in the case 
of v e — > u^. No day/night effect should be present, but a small seasonal effect might be 
expected (vacuum solution). 

The atmospheric neutrino anomaly should be confirmed, including the azimuthal depen- 
dence of the \xje suppression. Since it is mainly due to disappearance into sterile neutrinos, 
the muon deficit (rather than electron excess) should be confirmed, and the NC/CC ratio 
should not be enhanced. Experiments looking at upward going muons JO] should also see 



an effect, and current reactor experiments testing — v e oscillations p5]|[53| should continue 



to observe negative or marginal results. Long baseline accelerator experiments 
should not observe a strong v T appearance effect, but should observe maximal disappear- 
ance if the — uz oscillations lie in the experimentally accessible range. Unfortunately, the 
low Am 2 range (Am 2 < 10~ 3 eV 2 ) allowed by Superkamiokande |32 can not be fully tested 
by any of the currently planned long baseline experiments. 



However, ICARUS |5|| and MINOS |p5] , if sensitive enough, should observe energy in- 
dependent — v e oscillations, whose (small) rate should finally settle the value of sin 2 2# for 
the LSND result. 

~ v t oscillations should be observed in short baseline accelerator experiments within 
the next decade. If not by the current experiments NOMAD and CHORUS [J45| , then by 
future more sensitive experiments |57] with a slightly lower Am 2 threshold, — v T solution 
(b) would also be within the reach of the long baseline r appearance experiments |54|[56], 
and could possibly be confused with a signal at low Am 2 in the atmospheric neutrino region. 
Therefore both short and long baseline experiments are needed to unambigously confirm or 
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rule out this oscillation scenario. Since is of order eV or more (hierarchy), r neutrinos 
should make up a significant part of the missing dark matter. 

If — v T oscillations would be discovered with sizeable mixing, v T — disappearance 
could be verified in the long term by comparing the r production rates in short, medium, and 
long baseline experiments. Some corrections to the atmospheric neutrino spectrum might 
also be observable. 

Lastly, if this pseudo-Dirac neutrino solution would be confirmed and some of its param- 
eters measured more precisely, the complicated interplay between lepton number violating 
active-sterile neutrino oscillations, matter effects, and Big Bang nucleosynthesis (see ap- 
pendix) could be used to further constrain the remaining parameters. 



Appendix: Compatibility with Big Bang nucleosynthe- 
sis 

As is well known |58| , the ratio of deuterium to hydrogen (D /H) and the fraction of He in 
the universe (Y) is determined in the Big Bang model by the ratio of neutrons to protons 
at the time of the "weak freeze out" , i.e. the time when the reactions n + e + <-» p + v e 
and p + e~ <-> n + v e proceed too slowly compared to the expansion rate of the universe to 
keep n/p at its thermal equilibrium value. The effective number of light neutrino flavours 
N v enters as a contribution to the energy density, which in turn influences the expansion 
rate. From this effect, a typical upper limit derived from current Big Bang nucleosynthesis 
(BBN) models is N v < 3.1 @ 95% c.l. [5SJ with a central value of about 2.3, only marginally 



compatible with the standard model N v = 3. 

Active-sterile neutrino mixing can influence this rate in two ways: If active-sterile os- 
cillations occur during thermal equilibrium, i.e. significantly before weak freeze out, the 
sterile neutrinos will make an additional contribution to the energy density at freeze out, 
and therefore increase the expansion rate, effectively behaving like additional light neutrino 
generations. On the other hand, depletion of u e due to v e — > vf during freeze-out reduces 
the reaction rate for n/p interchange. In both cases, the predicted D/H and Y values will 
be larger than in the standard N v = 3 case, and therefore even less compatible with the 
observed values which yield the limit quoted above. 

The Am 2 of 10 -10 eV 2 or lower for vacuum v e — oscillations yields an oscillation time 
scale which is significantly longer than the freeze out time scale. The oscillation effect is thus 



negligible for practical purposes ||11|| . Both — and v T — oscillations approximately 



fall in the mass range 10 — 10 eV for which the BBN implications are discussed in ref. 



12| . Ignoring non standard model effects such as lepton number violation, each of these 
contribute 1 unit to N v if Am 2 is at the upper edge of this range, and slightly less than one 
unit if Am 2 is of order 10" 3 - 10" 4 Q, yielding N v ~ 4.6 - 5. 

However, matter effects induced by a significant initial lepton number asymmetry (10~ 3 
or larger) can almost fully suppress the — and v r — oscillations during thermal 
equilibrium |59" 1, therefore avoiding any contribution to the effective number of neutrinos. 



Moreover, such an asymmetry can be created by the (lepton number violating) active- 
sterile oscillations themselves ||60|| ) if small inter-generation oscillations (i.e. v T — ic) are 
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added. The quantitative treatment of the — u~[ case in ref. pT]| is qualitatively also 



A 4 "\l 

applicable to v T — oscillations. Interestingly, the parameters needed to create a lepton 
number asymmetry which suppresses — and v T — z/f oscillations until they become 
irrelevant for BBN largely overlap with the parameter space suggested by the pseudo-Dirac 
solution (table ||). For example, taking Am 2 j^ II2 = 1CT 3 eV 2 , Am 2 2 = 1 eV 2 (e = 3 x 1CT 3 
for case (1)) and adding v T — or v T — oscillations with sin 2 2# ~ 10~ 6 , which have 
a negligible effect on the phenomenology of the pseudo-Dirac model, almost full vx — Vx 
suppression (X = /i, r) is obtained for Am 2 3 > 3 eV 2 , which is just what is required from 
the hierarchy assumption. 

Finally, with a similar set of parameters, the "reprocessing" of the initial lepton number 
asymmetry into an effective excess of v e over i> e and a corresponding asymmetry in the 
n —>■ p and p — > n reaction rates during BBN can induce a reduction of the effective number 
of neutrinos by as much as 0.5 |62| . This yields the lower bound N v > 2.5, compatible with 
the current limit. Furthermore, the systematic error on the experimental input for Y and 
D/H is a matter of ongoing discussion fl63"l , which could contribute to loosening the BBN 
constraint. 

In conclusion, the apparent conflict between the BBN limit and 3 additional sterile neu- 
trinos can be evaded through the lepton number violating (i.e. non standard model) effects 
caused by active-sterile neutrino oscillations. 
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